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ABSTRACT 

Air  preheater  is  one  of  the  accessories  of  steam  generators,  that  works  as  a  heat  exchanger  that  utilizes  the  heat 
of  the  flue  gas  of  the  boiler,  heating  the  air  needed  for  the  combustion  of  the  fuel,  thereby  increasing  the  thermal 
efficiency  of  the  boiler.  It  makes  a  considerable  contribution  to  improve  the  overall  efficiency  of  fossil-fuel-fired  power 
plants. 

The  paper  presents  the  design  of  a  Recuperative  Tubular  type,  air-through-air  Preheater  with  the  staggered  tube 
arrangement,  designed  to  meet  specific  performance  requirements,  using  the  software  application  SOLIDWORKS. 
An  analytical  study  was  undertaken  and  various  heat  transfer  performance  parameters  like  outlet  temperature  of  the  air, 
inside  and  outside  heat  transfer  coefficients,  overall  heat  transfer  coefficient  and  rate  of  heat  transfer  involving 
conduction  and  convection  modes  of  heat  transfer  were  found.  The  results  were  tabulated  and  graphically  depicted. 
These  results  were  validated  using  CFD  analysis  at  100%  load  and  were  found  to  be  in  good  agreement  with  each  other, 
with  very  little  deviation.  These  heat  transfer  parameters  are  critical  in  the  designing  and  functioning  of  Air  Preheater  to 
get  the  desired  output.  Air  flow  phenomenon,  flue  gas  flow  phenomena,  temperature  distribution  and  pressure  drop  in  the 
Air  Preheater  were  discussed  in  the  CFD  analysis. 

KEYWORDS:  Air  Preheater,  Outlet  Temperature  of  Air,  Logarithmic  Mean  Temperature  Difference  LMTD, 
Nusselt  Number,  Prandtl  Number,  Reynolds  Number,  Overall  Heat  Transfer  Coefficient  &  Computational  Fluid 
Dynamics(CFD) 
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INTRODUCTION 

Air  Preheaters  are  designed  to  extract  the  hot  exhaust  energy  which,  when  discharged  unused  into  the 
atmosphere,  could  be  wasted.  Reduction  of  flue  gas  exit  temperature  helps  in  the  reduction  of  harmful  gases  up  to  a 
certain  extent.  So,  Air  Preheaters  contribute  to  controlling  environmental  pollution  and  global  warming. 
Air  Preheaters  help  in  reducing  coal  consumption,  which  helps  to  counter  the  rising  fuel  prices  and  make  a 
considerable  contribution  to  improve  the  thermal  efficiency  of  boilers  and  overall  efficiency  of  fossil  fuel  fired 
power  plants.  Air  Preheater  is  installed  between  the  economic  and  chimney  of  the  steam  power  plant  as  shown  in 
Figure  1. 


www.tjprc.ors 


SCOPUS  Indexed  Journal 


editor@tjprc.  org 


Original  Article 


1010 


J  Emeema,  MD  Khimran  Uddin  &  MD  Ashfan  Pasha 


Figure  1:  Schematic  Diagram  of  Plant  Layout 


Recuperative  Tubular  type  air  preheaters  have  an  enclosed  rectangular  box,  type  construction,  consisting  of 
several  straight  tubes  connecting  two  thick  end  plates  at  two  sides  of  the  box.  They  consist  of  straight  tube  bundles  which 
pass  through  the  outlet  ducting  of  the  boiler  and  open  at  each  end  outside  of  the  ducting.  Inside  the  ducting,  the  hot  furnace 
gases  pass  around  the  preheater  tubes,  transferring  heat  from  the  exhaust  gas  to  the  air  inside  the  tubes.  Ambient  air  is 
forced  by  a  fan  through  the  ducting  at  one  end  of  the  preheater  tubes  and  at  the  other  end,  the  heated  air  from  inside  of  the 
tubes  emerges  into  another  set  of  ducting,  which  carries  it  to  the  boiler  furnace  for  combustion. 


In  staggered  arrangement,  multiple  rows  of  tubes  are  arranged  in  planes  perpendicular  to  the  airflow  direction. 
Tubes  in  one  plane  are  spaced  at  a  pitch  of  “air  side  pitching”  (also  called  transverse  pitching).  Next  row  of  tubes  are 
placed  along  the  airflow  direction,  at  half  the  “gas  side  pitching”  (also  called  longitudinal  pitching)  but  staggered  to 
previous  row  of  tubes  by  half  of  airside  pitch.  Third  row  is  the  repeat  of  T'row  and  fourth  row  is  the  repeat  of  2ndrow  and 
so  on.  This  arrangement  is  repeated  to  get  the  total  number  of  tubes  per  pass. 


DESIGN  OF  AIR  PRE-HEATER 


The  Air  Preheater  designed  for  the  analysis  is  the  Recuperative  Tubular  type  with  the  staggered  tube  arrangement. 
It  consists  of  1488number  of  steel  tubes/pipes  arranged  in  32  rows,  with  47  tubes  in  1 sl  row  and  46  tubes  in2ndrow  and  so 
on.  Either  side  of  the  tubes  is  welded  into  the  plate.  It  is  designed  such  that  air  flows  through  the  tubes  and  flue  gas  flows 
over  the  tubes  heating  the  tube  surface.  The  design  of  Air  Preheater  is  done  in  SOLIDWORKS  which  is  a  comprehensive 
suite  of  computational  fluid  dynamics  software  for  modeling  fluid  flow  and  other  related  physical  phenomena. 

Design  Notations  and  Specifications  of  Air  Preheater 

Inside  diameter  of  tube  D  =  25mm 
Outside  diameter  of  tube  Dc)  =  28mm 
Length  of  tube  L  =  1.66m 
Width  of  tube  W  =  2.35m 


Impact  Factor  (JCC):  7.6197 
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Number  of  tubes  n  =  1488 

The  design  and  working  process  of  air  Preheater  is  shown  in  Figure  2. 


Figure  2:  Design  and  Working  Process  of  Air  Pre-Heater 


HEAT  TRANSFER  TERMINOLOGY 

Logarithmic  Mean  Temperature  Difference  LMTD 

The  LMTD  is  a  logarithmic  average  of  the  temperature  difference  between  the  hot  fluid  (flue  gas)  and  cold  fluid 
(air)  at  each  end  of  the  Air  Pre-heater.  The  larger  the  LMTD,  the  more  the  heat  transferred. 

LMTD  =  (AT  rAT2)  /  LN  x  (ATi/  AT2) 

ATi  =  Tgi  -  Tao  =  Temperature  difference  between  the  two  streams  at  end  1, 

AT2  =  Tg0  -  Tai=  Temperature  difference  between  the  two  streams  at  end  2. 

Tgi  =  Temperature  of  flue  gas  at  inlet  of  air  Preheater 

Tg0  =  Temperature  of  flue  gas  at  outlet  of  air  Preheater 

Tai  =  Temperature  of  air  at  inlet  of  air  Preheater 

Tao  =  Temperature  of  air  at  outlet  of  air  Preheater 

For  Cross-flow  heat  exchanger,  (LMTD)Crossflow=  (LMTD)Counterflow  x  F 


F  =  Correction  Factor  obtained  from  a  heat  transfer  data  book 
Dimensionless  Numbers  in  Forced  Convection 
Nu  =  Nusselt  Number 
Re  =  Reynolds  Number 
Pr  =  Prandtl  Number 


Overall  Heat  Transfer  Coefficient  Based  on  Outside  Surface 
U0  (with  Fouling  Factor) 


—  =  +  Rfo  +  ^  In  —  +  —  Rfi  +  —  x 

u0  ho  k  q  rj  rj  hj 


Rf0=  Flue  gas  side  or  outside  Fouling  Factor 
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Rfi=Air  side  or  inside  Fouling  Factor 

h„=  Flue  gas  side  or  outside  heat  transfer  coefficient 

hj=  Air  side  or  inside  heat  transfer  coefficient 

r0=  Outside  radius  of  the  tube 

r,  =  Inside  radius  of  tube 

k=  Thermal  conductivity  of  tube  material 

Overall  Heat  Transfer  Coefficient  Based  on  Outside  Surface  Area 
U0  (without  Fouling  Factor) 

^=±+£oinEo+£ox  i 
U0  h0  k  q  rj  hj 

Note:  Air  Preheater  is  designed  such  that  air  flows  inside  the  tubes  and  flue  gas  flow  outside  the  tubes.  Hence  the 
flow  parameters  are  denoted  as  ‘inside’  for  air  side  and  ‘outside’  for  flue  gas  side. 

Rate  of  Heat  Transfer  Q=  U0A0(LMTD)Crossflow 

A()  =  Outside  surface  area  of  the  tube  =  7iD0Ln 

ANALYTICAL  ANALYSIS  OF  AIR  PRE-HEATER 

Input  Parameters  of  Air  Preheater  At  100%  Load 


Inlet  temperature  of  flue  gas  T„j:  180°C 

Outlet  temperature  of  flue  gas  Tgo:130°C 
Inlet  temperature  of  air  Tai  30°C 

Volumetric  flow  rate  of  air  10.417  m3/s 

Volumetric  flow  rate  of  Flue  Gas  13.88  m3/s 
Tube  wall  temperature  Tw  120°C 


The  analysis  is  performed  at  three  different  loads  (40%,  100%,  and  120%)  on  the  air  preheater. 
Note:  The  values  of  input  Parameters  vary  for  different  loads. 

Output  Parameters  Found 

•  The  outlet  temperature  of  Air,  T.in 

•  Inside  Heat  transfer  coefficient  (air  side)  h; 

•  Outside  heat  transfer  coefficient  (flue  gas  side)h0 

•  Overall  heat  transfer  coefficient  based  on  outside  surface  area  U0. 

•  Rate  of  heat  transfer  Q 


Impact  Factor  (JCC):  7.6197 


SCOPUS  Indexed  Journal 


NAAS  Rating:  3.11 


Design  and  Analysis  of  Air  Preheater 


1013 


Heat  Transfer  Analysis  on  Flue  Gas  Side 

Effective  Area  of  passage  of  flue  gases 


A  =  W  x  L-  (No.  of  tubes  in  first  row  x  D0x  L) 


Velocity  of  Flue  Gas  V  = 


Volumetric  flow  rate  of  flue 
Effective  area  of  flue  gas 


gas 


Bulk  Temperature  (Tb)=  (Tgi+Tg0)/2 
Film  Temperature  (Tf)=  (Th+Tw)/2 

Properties  of  Flue  gas  taken  at  (Tf)  from  the  Heat  Transfer  Data  Book 
Density  (p) 

Kinematic  Viscosity  (v) 

Prandtl  Number  (Pr) 

Thermal  Conductivity  (k) 

Specific  heat  (Cpg) 

R  -VDo 

The  value  of  the  Reynolds  number  indicated  that  the  flow  on  the  flue  gas  side  is  Turbulent 

Convective  Heat  Transfer  Correlation  between  the  dimensionless  numbers  obtained  from  Heat  Transfer  Data 
Book  for  Turbulent  flow  in  forced  convection  is  as  follows. 

Nusselt  number  Nu  =  0.33  x  Re0  6  x  Pr0,33 


Nu  = 


hoxDo 

k 


hG  =  Flue  gas  side  heat  transfer  coefficient  or  outside  heat  transfer  coefficient  was  obtained. 
Heat  Transfer  Analysis  on  Air  Side 

An  energy  balance  between  hot  and  cold  fluids  gives: 


Q  =  mgx  Cpg(Tgo-  Tgi)  =  ma  x  Cpa  (Tao  -  Tal) 

Where  mg  =  mass  flow  rate  of  flue  gas  =  Volumetric  flow  rate  of  the  flue  gas  x  density  of  flue  gas 
ma  =  mass  flow  rate  of  air  =  Volumetric  flow  rate  of  air  x  density  of  air 
Tao  =  Temperature  of  air  at  outlet  is  obtained  from  the  above  equation 

T  +T  ■ 

Bulk  mean  temperature  of  air  Tb  =  a°2  al 

Properties  of  air  taken  at  (Tb)  from  Heat  Transfer  Data  Book 
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Density  (p) 

Kinematic  Viscosity  (v) 
Prandtl  Number  (Pr) 


Thermal  Conductivity  (k) 

Specific  heat  (Cpa) 

R  -VDi 

The  value  of  the  Reynolds  number  indicated  that  the  flow  on  the  air  side  is  Turbulent 


Convective  Heat  Transfer  Correlation  between  the  dimensionless  numbers  obtained  from  Heat  Transfer  Data 
Book  for  Turbulent  flow  in  forced  convection  based  on  the  value  of  the  Reynolds  number  and  L/D  ratio  is  as  follows. 

Nusselt  number  Nu  =  0.023  x  Re0  8  x  Pr0  4 


Nuxk 


h;  =  air  side  heat  transfer  coefficient  or  inside  heat  transfer  coefficient  was  obtained. 


Overall  Heat  Transfer  Coefficient  Based  on  Outside  Surface  Area 
U0  (with  Fouling  Factor) 


—  r~+  Rfo  +  -p-  In  —  +  —  Rfi  +  —  x 
u0  h0  k  !■[  !■[  rj  hj 


Rfo  =  0.002 
Rfi  =  0.0001 


Overall  Heat  Transfer  Coefficient  Based  on  Outside  Surface  Area 
U0( without  Fouling  Factor) 

^=±+b>lnEo+fox  1 
u0  h0  k  q  r-j  hj 

Rate  of  Heat  Transfer  Q  =  UD  A0(LMTD)Crossflow 


CALCULATED  VALUES  IN  TABULAR  FORM 


Table  1:  Flue  Gas  Side  Values  at  Different  Loads 


Load 

Volumetric  Flow  Rate  m  Vs 

Tgo  J 

Re 

Pr 

Nu 

h„ 

\\  7  in 2  K 

40% 

5.55 

150 

3442.442 

0.6815 

38.519 

48.19 

100% 

13.88 

130 

8795.3783 

0.683 

67.597 

83.528 

120% 

16.656 

135 

10490.6 

0.682 

75.187 

93.178 

Impact  Factor  (JCC):  7.6197 
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Table  2:  Air  Side  Values  at  Different  Loads 


Load 

Volumetric  Flow  Rate 

m3/s 

Tao1 

Re 

Pr 

Nu 

hi  W/m2  K 

40% 

5.55 

47.94 

8048.361 

0.698 

26.536 

29.82 

100% 

13.88 

89.31 

18798.46 

0.696 

52.24 

60.51 

120% 

16.656 

83.67 

22985.9 

0.6966 

61.38 

70.52 

Table  3:  Output  Parameters  at  Different  Loads 


Load 

Volumetric  Flow  Rate 

m'Vs 

T 

A  ao 

h„ 

W/m2  K 

hi 

W/m2  K 

u0 

W/m2  K 

Q 

kW 

40% 

5.55 

47.94 

48.19 

29.82 

17.12 

420.12 

100% 

13.88 

88.98 

83.53 

60.51 

30.66 

621.92 

120% 

16.656 

83.67 

93.17 

70.52 

34.78 

760.89 

CALCULATED  VALUES  IN  GRAPHICAL  FORM 


Load  vs  Tao 
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Graph  1:  Load  vs  Outlet  Temperature  of  Air  (Tao) 


Graph  2:  Load  vs  Overall  Heat  Transfer  Coefficient  (U0) 


Graph  3:  Overall  Heat  Transfer  Coefficient  (U0)  vs 
Outlet  Temperature  of  Air  (Tao)  for  Different  Loads 
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Re  vs  (air  side) 
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Graph  4:  Reynolds  Number  (Re)  vs  Inside  Heat  Transfer 
Coefficient  (hi)  for  Different  Loads 


Graph  5:  Reynolds  Number  (Re)  vs  Outside  Heat 
Transfer  Coefficient  (ho)  for  Different  Loads 


Graph  6:  Reynolds  Number  (Re)  on  Air  Side  vs  Overall  Heat 
Transfer  Coefficient  (Uo)  with  and  without 
Fouling  Factor  at  Different  Loads 

COMPUTATIONAL  FLUID  DYNAMICS  (CFD) 

ANALYSIS  OF  AIR  PRE-HEATER 


A  CFD  is  a  method  for  solving  complex  fluid  flow  and  heat  transfer  problems  on  a  computer.  CFD  allows  the 
study  of  problems  that  are  too  difficult  to  solve  using  classical  techniques.  The  flow  inside  the  ESP  is  complex  and  this  can 
be  analyzed  using  the  CFD  tool,  which  provides  an  insight  into  the  complex  flow  behavior. 

The  CFD  Analysis  is  performed  at  100%  load  to  find  the  following  parameters: 

•  Outlet  temperature  of  air  Tao 
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•  Outlet  temperature  of  Flue  Gas  T  go 

•  Overall  heat  transfer  coefficient  UQ 

•  Temperature  distribution  along  the  length  of  the  tube 

•  Pressure  distribution  along  the  length  of  the  tube 

Inlet  Boundary  Conditions  for  CFD  Analysis  is  Shown  in  Figure  3 


Figure  3:  Inlet  Boundary  Conditions  for  CFD  Analysis 


Input  parameters  for  CFD  analysis  (same  values  used  in  analytical  analysis) 

Inlet  temperature  of  Air  Tai  =  30°  C 

Inlet  temperature  of  flue  gas  Tgi  =  180°  C 

Volumetric  flow  rate  of  air  =  10.417  m3/s 

Volumetric  flow  rate  of  Flue  Gas  =  13.88  m3/s 


RESULTS  FROM  CFD  ANALYSIS 

Temperature  Variation  of  Air  from  Inlet 
to  Outlet  is  shown  in  Figure  4 


Figure  4:  Top  View  of  Air  Pre-Heater  with  Inlet  and  Outlet 
flow  Variation  of  Temperature  of  Air 


Blue  color  indicates  inlet  temperature  of  ambient  air 


Green  color  indicates  outlet  temperature  of  hot  air 
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Temperature  Variation  of  Flue  Gas  from  Inlet 
to  Outlet  is  shown  in  Figure  5 


Figure  5:  Front  View  of  Air  Pre-Heater  with  Inlet  and 
Outlet  flow  Variation  of  Flue  Gas  Temperature 

Velocity  Variation  of  Flue  Gas  from  Inlet 
to  Outlet  is  shown  in  Figure  6 


Figure  6:  Top  View  of  Air  Pre-Heater  with  Inlet  and 
Outlet  Flow  Variation  of  Velocity  of  Air 


Pressure  Variation  of  Air  from  Inlet 
to  Outlet  is  shown  in  Figure  7 


Figure  7:  Top  View  of  Air  Pre-Heater  with  Inlet  and 
Outlet  Flow  Variation  of  Relative  Pressure  of  Air 


Impact  Factor  (JCC):  7.6197 
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Graph  7:  Temperature  Distribution  of  Air  and 
Flue  Gas  across  the  Length  of  Air  Pre-Heater 


Graph  8:  Pressure  Distribution  across  the  Length  of 
Air  Pre-Heateris  shown  in  Graph  8 


Length  (m)  vs  Pressure  (Pa) 


101325.01  - 

101325.01 

j; 

101325 

0  0.5  1  1.5  2  2.5 

Length  [m] 


Graph  8:  Pressure  Distribution  across  the 
Length  of  Air  Pre-Heater 


Graph  8  Legend 

Black  line  indicates  the  flue  gas  pressure 
Gray  line  indicates  the  air  pressure 

COMPARISON  BETWEEN  ANALYTICAL  AND 
CFD  ANALYSIS  AT  100%  LOAD 

The  output  Heat  Transfer  parameters  obtained  from  Analytical  Analysis  and  CFD  Analysis  are  Tabulated  in 
Table4  for  comparison. 


Table  4:  Values  of  Output  Heat  Transfer  Parameters 
from  Analytical  &  CFD  Analysis 


Parameter 

Analytical  Analysis 

Cfdanalysis 

Outlet  temperature  of  air  Tao  in  □ 

88.98 

88.07 

Outlet  temperature  of  flue  gas  TB0  in  □ 

130 

127.33 

w 

Overall  heat  transfer  coefficient  U0(^^) 

30.66 

31.82 

Rate  of  heat  transfer  Q  in  kW 

621.92 

645.45 

To  find  the  overall  heat  transfer  coefficient  U0in  CFD  analysis,  the  outlet  Temperatures  of  Air  and  Flue  gas 
obtained  from  CFD  result  are  used  to  calculate  the  values  of  inside  and  outside  heat  transfer  coefficients  h,  and  h0  and 
overall  heat  transfer  coefficient  U0. 
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RESULTS  AND  DISCUSSIONS 

Analytical  Analysis 

The  following  observations  can  be  made  from  graphs. 

•  Graph  3:  UG  increased  with  load.  Tao  is  maximum  at  100%  load 

•  Graphs  4  &  5:  hj  and  hQ  increase  with  Reynolds  number  Re 

•  Graph  6:  UD  reduces  due  to  Fouling  by  6.5  % 

•  From  Table  3,  h,<  h0  which  indicates  that  the  air  side  heat  transfer  coefficient  offers  more  thermal  resistance  and 
controls  the  rate  of  heat  transfer.  Therefore,  some  enhancement  technique  like  fins  can  be  employed  on  the  air 
side  to  compensate  for  low  values  of  h,. 

Comparison  of  Analytical  &CFD  Analysis  At  100%  Load 

•  The  outlet  temperature  of  air  Tao  increased  by  58.98°C  in  Analytical  analysis  and  by  58.07°  C  in  CFD  analysis. 

•  The  Overall  heat  transfer  coefficient  UG  from  CFD  analysis  has  a  deviation  of  3.7  %  of  Analytical  analysis 

•  The  outlet  temperature  of  air  Taofrom  CFD  analysis  has  a  deviation  of  1.03  %  of  Analytical  analysis 

•  Outlet  temperature  of  flue  gas  Tgo  from  CFD  analysis  has  a  deviation  of  2.05  %  of  Analytical  analysis 

•  The  rate  of  heat  transfer  Q  from  CFD  analysis  has  a  deviation  of  3.78  %  of  Analytical  analysis 
Temperature  Distribution  Graph  in  CFD  Analysis 

From  graph  7,  it  can  be  observed  that: 

•  The  temperature  of  ambient  Air  has  increased  by  57.98°  C,  from  303.09  K  (30.09  °C)  to  361.07  K  (88.07  °C) 

•  The  temperature  of  flue  gas  has  decreased  by  52.73°  C,  from  453.06  K  (180.06  °C)  to  400.33  K  (127.33  °C) 

Pressure  Distribution  Graph  in  CFD  Analysis 

Graph  8  indicates  that  the  pressure  of  flue  gas  along  the  length  of  Air  Preheater  has  dropped  by  0.01  Pa,  from 
101325.01  Pa  to  101325  Pa.  Hence,  it  can  be  concluded  that  the  pressure  drop  along  the  Air  Preheater  is  very  less. 

CONCLUSIONS 

From  the  design  and  analysis  of  Air  Preheater,  it  can  be  concluded  that  the  heat  transfer  performance  parameters 
like  outlet  temperature  of  air  Tao,  outlet  temperature  of  flue  gas  Tgo,  overall  heat  transfer  coefficient  based  on  outside 
surface  area  UQ  and  rate  of  heat  transfer  Q  obtained  from  analytical  analysis  are  in  good  agreement  with  the  corresponding 
values  obtained  from  CFD  analysis,  with  very  little  deviation.  Hence  it  can  be  concluded  that  these  design  parameters  of  air 
Preheater  met  with  the  performance  requirements  as  intended.  It  can  also  be  concluded  that  the  designed  air 
Preheaterserves  the  purpose  of  extracting  heat  from  the  flue  gas  of  boilers,  which  can  be  used  to  preheat  the  air  required  for 
the  combustion  of  the  fuel  in  the  boilers,  thereby  increasing  thermal  efficiency  of  the  boiler,  ultimately  making  a 
considerable  contribution  to  improve  the  overall  efficiency  of  fossil-fuel-fired  power  plants. 
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